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ABSTRACT 

The iron Ka line commonly observed in the X-ray spectrum of both stellar- 
mass and supermassive black hole candidates originates from X-ray fluorescence 
of the inner accretion disk. Accordingly, it can be used to map the spacetime 
geometry around these objects. In this paper, we extend previous work using 
the iron Ka line to test the Kerr black hole hypothesis. We adopt the Cardoso- 
Pani-Rico parametrization and we test the possibility of constraining possible 
deviations from the Kerr solution that can be obtained from observations across 
the range of black hole spins and inclination angles. We confirm previous claims 
that the iron Ka line is potentially a quite powerful probe for testing the Kerr 
metric given sufficiently high quality data and with systematics under control, 
especially in the case of fast-rotating black holes and high inclination angles 
since both conditions serve to maximize relativistic effects. We hnd that some 
geometric perturbations from Kerr geometry manifest more strongly in the iron 
line prohle than others. While the perturbation parameter e\ can be well con¬ 
strained by the iron line prohle, an orthogonal data set is necessary to constrain 
departures from Kerr geometry in €3. 

Subject headings: accretion, accretion disks — black hole physics — gravitation 


1. Introduction 

When a star exhausts its nuclear fuel, the core shrinks to hnd a new equilibrium con- 
hguration. In the case of a sufficiently massive star, there is no known mechanism capable 
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of providing the necessary pressure support to halt the process, and the body undergoes 
complete collapse. In general relativity and under a set of reasonable assumptions, the final 


product of the collapse is a black hole (BH) (Baiotti et al. 2005 Baiotti & Rezzolla 2006 


Giacomazzo et al. 2011; Joshi & Malafarina 2011[ ). In 4 dimensions, the only uncharged 
BH solution of the vacuum Einstein’s equations is the Kerr metric. The spacetime around 
astrophysical BHs should be well described by the Kerr geometry, since the presence of the 
accretion disk is usually completely negligible, given that the disk mass is many orders of 


magnitude smaller than that of the BH (Bambi et al. 2014) 


The objects termed BH candidates are dark and compact and are naturally interpreted 
as the Kerr BHs of general relativity. Many of these sources have robust mass measurements, 
obtained dynamically (e.g., Ozel et al. 2010, and references therein). When a compact 


object in an X-ray binary is more massive than ^ 3 Mq, the object is classified as a BH 


candidate, since it exceeds the maximum mass for a neutron star (Rhoades & Ruffini 1974). 


With a similar spirit, the supermassive bodies at the center of every normal galaxy are BH 
candidates because they are too massive, compact, and too old to be clusters of neutron 
star^ The absence of thermal radiation emitted from the surface of BH candidates may 


also be interpreted as evidence for the presence of a horizon ( 

Narayan & McGlintock 

2008 

Broderick et al. 

2009 

i but see 

Abramowicz et al. 

2 (1(12 

Bambi 

2012 b 

). Eventually, the Kerr 


paradigm relies on the validity of general relativity in strong gravitational helds. However, 
deviations from the Kerr solution may be expected from a number of theoretical arguments, 
see, e.g. Dvali & Gomez ( 2011| 2013); Giddings (2014). The experimental conhrmation of 
the Kerr nature of BH candidates is thus an important issue in the fields of astrophysics. 


gravity, and high energy physics (Bambi 2011, 2013d). 


The spacetime geometry around BH candidates can be potentially tested by studying 
the properties of the electromagnetic radiation emitted by the gas in the accretion disk. Such 
radiation is inevitably affected by relativistic effects and therefore carries information about 
the metric. Accurate measurements of the thermal spectrum of a thin disk (Zhang et al. 


1997 McGlintock et al.||2014 ) and of the iron Ka line ( Fabian et al.||1989 ; Reynold^|2014 ) are 


today the principal techniques for providing information about the nature of BH candidates 
and they can be used to test the Kerr metric (Bambi fc Barausse||2011 Bambi||2012a, 2013c 


2012 dle 

2013b 

2014b 

Johannsen & Psaltis||2013 

Torres||2002 

Lu & Torres||2003 


the usual problem with such a test is often strong degeneracy among system parameters, in 
particular between the parameter measuring deviation from Kerr geometry, and the spin and 
the inclination angle. The thermal spectrum of thin disks has a simple shape and therefore 


Mhe cluster lifetime due to evaporation and collisions is less than the age of these systems (Maoz 


19981. 
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it is fundamentally more challenging to break such a degeneracy (Kong et ah 2014 Bambi 


2014a). The iron line profile is potentially more powerful. With available X-ray data, it 


has already been established that BH candidates cannot be some types of boson stars and 


of traversable wormholes (Bambi 2013a Bambi & Malafarina 2013). However, it isn’t clear 


whether present data are sufficiently precise to test the Kerr metric; instead, this is a target 


for the next generation of X-ray satellites (Jiang et al. 2015). 


The aim of this paper is to extend previous studies which test the Kerr metric using of 
the iron Ka line. In what follows, we employ the simplest version of the Cardoso-Pani-Rico 


parametrization and we consider two types of deformation in their schema, Cg and Cg (Cardoso 
et al. 2014). We discuss constraining eg and Cg using iron lines for different values of the BH 


spin and for a range of disk inclination angles with respect to our line of sight. We find that 
the best source with which to test the Kerr metric has high spin and a large inclination. 
Indeed, for a fast-rotating Kerr BH, the inner edge of the disk extends very close to the 
compact object, increasing the iron line’s ability to probe the region in which deviations from 
the Kerr geometry are most likely to be pronounced. A large inclination angle can likewise 
maximize the light bending effects, and accordingly amplify any differences resulting from 
alternate spacetimes. Some modes of deformation are surely more readily constrained than 
others. The deformation parameter e\ principally alters the position of the inner edge of the 
disk, whose effects on the iron line profile is generally more prominent than a non-vanishing 
Cg, whose main effect is to upon photon trajectories. The best candidates source to test for 
these variations would be a system like GRS1915-I-105, which is quite bright and harbors 
what appears to be a fast-rotating Kerr BH at a relatively large inclination angle (a* > 0.98 
McClintock et al. 2006| Reid et al. 2014). 


and i = 60° 


Throughout the Paper, we employ units in which = c = 1 and the convention of a 
metric with signature (—h ++). 


2. Testing the Kerr paradigm 


Gas in the accretion disk and the photons emitted by the disk can be treated as test- 
particles and are assumed to follow the spacetime geodesics; this is always the case in a 
metric theory of gravity but is not so in more general frameworks. The adopted approach is 
accordingly only sensitive to the geodesics of the metric and it cannot completely test the 
Einstein equations, in the sense that it is fundamentally impossible to use such a technique to 
distinguish a Kerr BH of general relativity from a Kerr BH in an alternative metric theory 


of gravity (Psaltis et al. 2008). Bearing this in mind, the usual strategy is to employ an 


approach similar to the Parametrized Post-Newtonian (PPN) formalism (Will 2006), which 
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has been used to test the Schwarzschild solution in the weak held limit within the Solar 
System since the 1960s. In the PPN case, one writes the most general static and spherically 
symmetric solution in a weak gravitational held, by using an expansion in M/r. In isotropic 
coordinates, the line element turns out to be 


ds^ = - 1- 


2M 2M^ 
- +/ 3 ^ + 


dt‘ 


+ {dx^ + dy^ + dz^) , 


( 1 ) 


where fd and 7 parametrize our ignorance, while the term 2M/r in is necessary to have the 
correct Newtonian limilQ In general relativity, we know that the only static and spherically 
symmetric vacuum solution is the Schwarzschild metric, and cast in the form above it has 
/3 = 7 = 1 . Since we want to test the Schwarzschild metric, we leave fd and 7 free and 
we aim to constrain the two parameters with observations. Today’s best measurements 


are (Williams et al. 2004 Bertotti et ah 2003) 


I/? - 1| < 2.3 ■ 10-^ |7-1|<2.3-10 


-5 


(3) 


To test the Kerr metric, we can use a similar approach, namely we consider a more 
general solution that includes the Kerr metric as special case. Such a metric will be described 
by a set of “deformation parameters”, whose value will be determined by observations and 
the Kerr geometry is recovered when all the deformation parameters vanish. However, now 
we want to test the Kerr metric near the BH, so we cannot use an expansion in M/r any 
longer. Spacetime is axisymmetric, not spherically symmetric, and if we try to alter by hand 
the Kerr metric it is difficult to avoid creating a spacetime with pathological features like 
naked singularities and regions with closed time-like curves. For the time being, there is not 
yet a standard formalism for testing the Kerr metric, even if some proposals are present in 


the literature (Johannsen & Psaltis 2011 Cardoso et al. 2014 Rezzolla & Zhidenko 2014). 


In what follows, we adopt the Cardoso-Pani-Rico parametrization (Cardoso et al. 2014). In 


^The PPN approach is traditionally formulated in isotropic coordinates and the line element of the most 
general static and spherically symmetric distribution of matter is given by Eq. Q. When cast in the more 
familiar Schwarzschild coordinates, the line element becomes 


ds^ = -[l- 


2M , 2M2 \ , 2 2M > , 2 ,^2 2 . 2 ^ , ,2 

-h (/3 - 7 ) —^ + ...] dt^ + [1+ 7 -h ... dr^ + r^de^ + sm^ Odcj) . (2) 
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Boyer-Lindquist coordinates, the line element reads 


, , I 2Mr\ , , . S(l + r) , . ^ . 

ds — — ( 1 -zq— ] (l + h ) (it H— -- —-—:—^^dv + YidO 

S y ^ ^ A + sin^ 9 


— 2 asin^ 6 

+ sin20|s + a2 sin^ 0 


y(i+ h‘) (1 + h’-) - (1 - (1 + h<) 


dtdcf) 


2a/( 1 + h*) (1 + h^) ~ -(l + ^*) 


(4) 


where a = J/M is the BH spin parameter, J is the BH spin angular momentum, S = 


+ cos^ 0, A = — 2Mr + a^, and 


= 


H-OO 

E 

k=0 


-2k 


+ 6 ' 


Mr 


2k+l' 


M2 

S 


^ , hJ- = 


+ 00 

E 

fc =0 


'-2k 


+ e 


Mr 


2k+l' 


M\ 

ifl ■ ® 


There are two infinite sets of deformation parameters, {e^} and {e^} (at increasingly high 
order). The line element in (Q reduces to the Johannsen-Psaltis one for h* = (Johannsen 
& Psaltis 2011), and to the Kerr line element when = 0. In the Johannsen- 


Psaltis background, Cq = 0 gives an asymptotically fiat spacetime, while ei and 62 must be 


small to meet the Solar System constraints (Cardoso et ah 2014). 63 is the first unbounded 


deformation parameter and there are no qualitative differences between 63 and higher order 
terms, in the sense that all these parameters deform the spacetime geometry in a similar way 
and leave the same observational signature in the spectrum ( Bambi||2012c ). In what follows, 
we only consider the possibility of non-vanishing Cg and Cg. These deformation parameters 
introduce deviations from the Kerr geodesic motion, which can be seen as a departure from 
the predictions of Einstein’s theory of gravity on how matter distorts the nearby spacetime. 


3. Simulations 


The broad iron Ka line observed in the X-ray spectrum of BH candidates of all masses 
originates when hard X-rays produced in a corona illuminate the accretion disk, fluorescing 


the line emission (Reynolds 2014). Geometrically thin, optically thick accretion disks can 


be described by the Novikov-Thorne model (Novikov fc ThornulllOTS; Page fc Thorne||1974), 


in which the gas particles follow geodesic circular orbits in the equatorial plane and the 
inner edge of the disk extends inwards to the innermost stable circular orbit, an assumption 
supported empirically by the stability of the measured inner edge of the disk across time and 
over an order of magnitude of luminosity variation (e.g., Steiner et al.||2010). Calculating the 


iron line profile for a non-Kerr background is straightforward and discussed elsewhere (e.g., 

Johannsen & Psaltis 2013). The effect of the deformation 


Bambi 

2012 dle, 

2013b, 

2014b 
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Fig. 1.— Impact of the deformation parameters Cg and Cg on the prohle of the iron line. 
The spin parameter is a* = 0.5 in the top panels and a* = 0.95 in the bottom panels. The 
viewing angle is always i = 45°. The left panels show the change which results from varying 
eg while keeping eg = 0. In the right panels, eg = 0 and we instead vary eg. 
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parameters e| and €3 on the iron line profile is shown in Fig. A non-vanishing deformation 
parameter mainly affects the extension of the low energy tail of the iron line prohle (as it 
changes the inner edge of the disk) as well as the photon redshift at 5-20 gravitational radii 
(which affects the shape of the iron line at high energies). 

In this section, we compnte the iron line assnming a Kerr backgronnd while estimating 
plausible constraints on e| and €3 to be obtained from observations. In doing so, we simulate 
the effect of observing our model in alternate backgrounds to assess how readily features can 
be distinguished at a given level of conhdence. In this work, we use the approach employed 

. We adopt a primed notation (e.g., a(, i') for the input model from 
simulated data are derived. Unprimed terms are used to refer to the model employed by the 
spectral htting routine. 


m 


Jiang et ah (2015 


Our simulation consists of hrst computing the photon flux of the iron line prohle over a 
grid of energy bins with resolution AE = 100 eV from 1-9 keV. For the sake of simplicity, 
the emissivity index for both the simulation and in the model htting is always assumed to be 
3; namely the intensity prohle scales as 1/r^. Such a value correspond to the Newtonian limit 
for lamppost coronal geometry (e.g., Fabian et ah 2012). A more exact model consisting of 


a radially evolving emissivity index - g - may be useful when treating real data (see, e.g., 
Dauser et al.||2013 ); however, in this exploratory work we adopt a widely-used simplihcation 
g = 3 in both generating our simulated data, and in its analysis. In our simulations, the 
iron line prohle is added to a power-law continuum. The latter is normalized to include 100 
times the number of iron line photons when integrated over the energy range 1-9 keV, and 
is generated using a photon index F' = 2. Depending on the line shape, this corresponds to 
a line equivalent width of EW ~ 370 — 440 eV. In the analysis below, the ratio between the 
continuum and the iron line photon hux, K, as well as the index of the continuum, F, are 
left as free parameters to be determined by the ht. 


We employ the standard analysis technique common throughout X-ray astronomy: our 
simulated spectra include Poisson noise, are then binned to achieve a threshold of counts, 
and htted employing analysis. The Poisson noise is trivial to include; in each bin it is 
a randomly generated realization using the “true” model prediction for the photon count. 
The simulated (noisy) photon count number, is what we employ as the data for our analysis. 
The level at which the noise manifests depends on the total number of spectral counts; in 
essence the total counts therefore set the scale of the noise in the spectrum. In this paper, 
we assume that the number of iron-line photons is 77 = lO'^, corresponding approximately 
to the line signal expected in a next generation high-quality observation of a supermassive 
black hole in an active galactic nucleus (AGN). To treat the simulation as if it were data, 
our analysis employs standard techniques, consisting of grouping the data to a minimum of 









20 counts per bin. Our results are not sensitive to this particular choice of threshold. 


Figs. [2]j^ show our results. We begin by considering the two deviation parameters 
separately, holding one fixed to zero, and allowing the other to vary. In Fig. the reference 
model is a Kerr BH (eg^ = Cg' = 0) with spin parameter a'^ = 0.5 and observed from a viewing 
angle F = 20° (left panels) and 70° (right panels). It is compared with a set of models in 
which the spacetime geometry is described by the Cardoso-Pani-Rico metric with only one 
non-vanishing deformation parameter, namely e\ in the top panels and Cg in the bottom 
panels. Here i, K, and F are always £t parameters and is minimized by varying them. 
These plots show reference contours of ~ Xmm = 1) 4, 9, 16, and 25. It is evident 

that Cg is more difficult to constrain from the line profile alone. A similar conclusion is found 
by Bambi (2014a). It is evident that iron lines are more capable of constraining deformation 
parameters when the inclination of the reference model i' is large, since this maximizes light 
bending effects. 


In Fig. the reference model is a fast-rotating Kerr BH with spin parameter a( = 0.95. 
As in Fig. the left panels are for an inclination angle i' = 20° and the right panels are 
for i' = 70°, the top panels for a non-vanishing deformation parameter e\ and the bottom 
panels for a non-vanishing deformation parameter Cg. In the right bottom panel, we only 
show Ay^ = 9, 16, and 25 for clarity, but we note that y^ has an almost degenerate valley 
extending to to large positive values of Cg. i is still free and y^ is minimized over i. We hnd 
that rapidly-rotating BHs make better targets for testing the Kerr metric, enabling the iron 
line analysis to produce more stringent constraints. Here also, a larger inclination angle is 
helpful, and again eg proves much more challenging to measure than eg. Even in the case 
of a fast-rotating BH with an almost edge-on disk, the iron line cannot really distinguish a 
Kerr spacetime from a background with a very large eg. A longer exposure time or a larger 
effective area of the detector, namely a higher number of photons A, would surely provide 
stronger constraints on eg, as y^ is approximately proportional to N. However, a higher 
photon count would not be very helpful for eg because the valley of y^ is almost degenerate: 
the area of the allowed region in Fig. would be smaller, but still large positive values of eg 
would be permitted. 


Fig. [^illustrates how the photon count number in the line A, which enters in the Poisson 
noise, affects the £t constraints. The values of the parameters are the same as in Fig. [^ but 
here we have A = 10^ photons in the iron line {K' = 100 as before). The constraints are 
markedly worse. A = 10^ photons corresponds approximately to the line signal available in 
a present-day high-quality observation of a supermassive black hole, while A = 10^ photons 
can be expected using the next generation of X-ray satellites. For instance, in the case of a 
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Fig. 2.— Top panels: contours with N = 10^ from the comparison of the iron line 

prohle of a Kerr BH simulated using an input spin parameter a'^ = 0.5 and an inclination of 
i' = 20° (left panel) and i' = 70° (right panel) versus a set of Cardoso-Pani-Rico BHs with 
spin parameters a* and a non-vanishing deformation parameter eg. Bottom panels: as in the 
top panels for a non-vanishing deformation parameter Cg. 
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Fig. 3.— Top panels: contours with N = 10^ from the comparison of the iron line 

prohle of a Kerr BH with an input spin parameter a'^ = 0.95 and an inclination of i' = 20° 
(left panel) and i' = 70° (right panel) versus a set of Cardoso-Pani-Rico BHs with spin 
parameters a* and a non-vanishing deformation parameter e\. Bottom panels: as in the top 
panels for a non-vanishing deformation parameter eg. 
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Fig. 4.— As in Fig. [^for A^ = 10^. 
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Fig. 5.— Impact of the power-law continuum. The reference model is a Kerr BH with 
a'^ = 0.95 and i' = 70° and the photon number in the iron line is = 10^. In the left panel, 
the simulation does not include the power-law component. The right panel, which is the 
same plot as the top right panel in Fig. we have both iron line and power-law component. 
The latter is normalized to include 100 times the number of iron line photons when integrated 
over the energy range 1-9 keV, and is generated using a photon index F' = 2. 
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Kerr BH with a'^ = 0.95 and i' = 70°, the contour Ax^ = 4 (9) provides the constraints 

-1 < 4 3 (-2 4 4 4 4.5) for fV = 10*. 

-0.5 <e‘,< 0.5 (-1 < 4 < 3) tor iV = 10“ . (6) 

The deformation parameter is unconstrained with both N = 10^ and N = 10^. However, 
for a hxed the uncertainty on the spin parameter is different. The contour = 9 yields 
a spin uncertainty Aa* ~ 0.3 for N = 10^, and Aa* ~ 0.05 for N = 10^. 

It is instructive to consider the effect of the power-law continuum on the constraining- 
power of the iron line signal. This is illustrated in Fig. [^for the deformation parameter e|. 
The reference model is a Kerr BH with a'^ = 0.95 and i' = 70°. The left panel employs the 
same model as in Fig. but with just the line, and no power-law component. The right panel 
is the same plot as the top right panel in Fig. The continuum’s presence adds additional 
photon noise to the line’s signal, which unsurprisingly weakens the resulting constraints, for 
a given line strength. The contour Ay^ = 25 would give the constraint 

0.9 < a, <1.1, -1< 6*3 <1.5 (7) 

without power-law component. The constraints become 

0.3 < a, <1.2, -l<e*3<4 (8) 


if we include the continuum. 

Lastly, Fig. shows the case in which both Cg and 63 are allowed to be non-vanishing 
at the same time. The reference model here is a fast-rotating Kerr BH with spin parameter 
a'^ = 0.95 and viewing angle i' = 70°. According to the previous plots, this is a fruitful 
source for testing the Kerr metric. Here, spin, inclination, K, and F are optimized in the £t. 
Strikingly, the measurements of 63 and 63 are essentially uncorrelated. 

The constraints obtained in Figs. [2]j^can be qualitatively understood as follows. When 
the spin parameter vanishes, gtrp = 0 , and the position of the inner edge of the disk is only 
determined by 63 . 63 affects instead the photon propagation from the disk to the distant 
observer. Since the position of the inner edge of the disk has a strong impact on the iron 
line prohle, 63 is much easier to constrain than 63 . For a non-vanishing spin, the picture is 
more complicated, but the general explanation is still valid and the influence of 63 on the 
inner edge of the disk is only moderate. 
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4. Discussion and conclusions 

An accurate measurement of the iron Ka iron line commonly observed in the X-ray 
spectrum of both stellar-mass and supermassive BH candidates is potentially a powerful tool 
for probing the spacetime geometry around these objects and test the Kerr BH paradigm. 
In this paper, we have extended previous studies by exploring how high quality X-ray data 
can test the nature of BH candidates for different BH spins, different viewing angles, and 
two independent types of deviations from Kerr. Assuming the case of a detection of = 10"^ 
photons in the iron line and Nc = 10® photons in power-law continuum, we have performed 
a number of simulations, whose constraints are reported in Figs. |2]j^ Our results can be 
summarized as follows: 

1. The best target sources to test the Kerr metric are BH candidates with high spin pa¬ 
rameter and observed from a large inclination angle. This is perfectly understandable, 
as both the situations enhance the magnitude of relativistic effects. 

2 . Observations can actually constrain the deformation parameter e\, while constraints 
on 63 seems to be out of reach with the sole use of iron line data. To measure 63 it 
is presumably necessary to combine the iron line measurement with a technique only 
sensitive to the spin parameter, and an example could be the determination of the 
spin through the jet power (e.g., as suggested by 

, if such a relationship is conhrmed. Since affects the photon propagation, it 
is also possible that reverberation data, which have a timing information, can help to 
constrain this parameter. 

3. While we have not yet explored the whole parameter space, our results suggest that 
the measurements of eg and of eg are not very correlated, namely it is possible to 
distinguish the effects of the two classes of deformation on the iron line prohle. 

In this exploratory work, we have not considered the effects of varying other key parameters 
related to Fe line modeling (emissivity index q, ionization parameter (, etc.). Moreover, we 
have made the simplifying assumption that systematic effects which may otherwise obscure 
our ability to assess these subtle deviations are under full control. In future work, it will 
be important to employ a more rigorous treatment of full reflection modeling, to assess 
systematic effects, and ultimately to perform a full and detailed analysis on the best available 
data of iron lines produced in strong gravity, in order to constrain any deviation terms 
directly. 
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Fig. 6.— Ax^ contours with N = 10^ from the comparison of the iron line prohle of a Kerr 
BH with a spin parameter a'^ = 0.95 and an inclination of i' = 70° versus a set of Cardoso- 
Pani-Rico BHs with non-vanishing deformation parameters e| and e^. See the text for more 
details. 




